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EXECUTIVE SUMMARY 


Clusters of metal atoms are of importance as models of substrates 
in studies of Chemisorption on metal surfaces and of catalysis of chemical 
reactions; they are also of increasing interest as novel molecular species 
in their own right that, surrounded by certain ligands, are possibly a whole 
new class of catalysts. Moreover, clusters of metal atoms are useful as a 
prototype atomic-level host systems fo’" the development and testing of models 
of interatomic forces. These interatomic forces are a basic ingredient to 
the much more complex atomic model "ng of the mechanical properties of metals 
that are useU to address import "c materials problems such as the effect of 
hydrogen on crack propagation. 

The approach we have taken to the study of metal atom clusters is 
based on the application of ab initio quantum mechanical approaches. This 
interim report discusses our current research effort in the application of 
these methods. Because these large "molecular" systems pose special practical 
computational problems in the application of the quantum mechanical methods, 
there is a special need to find simplifying techniques that do not compromise 
the reliability of the calculations. Our current research is therefore 
directed towards various aspects of the implementation of the Effective Core 
Potential (ECP) technique for the removal of the metal atom core electrons 
from the calculations. This interim report discusses our recent progress in 
this area. 

During this reporting period we have developed a reformulation of 
the ECP rrolecular integral formulae to bring out the shell characteristics 
in the formulae. This is an important first step towards fulfilling the 
increasing need, especially for the metal atom clusters, to speed up the 
computation of the ECP integrals. Also during this reporting period, we have 
completed and published work on tlie relationships among the derivatives of 
the molecular integrals that extends some of the techniques pioneered by 
Komornicki for the calculation of the gradients of the electronic energy. 
Finially, during this reporting period, we have discovered a new formulation 
of the ECP approach that quite naturally unifies the various state-of-the- 
art "shape- and Hamiltonian-consistent" techniques. 


QUANTUM CHEMICAL CALCULATION OF THE EQUILIBRIUM 
STRUCTURES OF SMALL METAL ATOM CLUSTERS 

by 

Luis R. Kahn 
INTRODUCTION 


The goal of this research is to gain quantitative fundamental 
understanding about the properties of clusters of metal atoms. The in- 
creasing interest in these novel "molecular" species stems from their use 
to model substrates in the study of chemisorption, and from the prospects 
that such studies will yield the long-sought understanding of the mecha- 
nisms of heterogeneous catalysis of chemical reactions. Moreover, these 
metal atom clusters surrounded by carbonyl groups, are indeed being thought 
of as new types of catalysts in their own right. In every case, however, 
there is still very little known yet about these interesting new molecular 
species. Since the atomic level properties of these clusters are not yet 
readily available from experiment, the most reliable approach to their 
determination lies in the development and application of predictive ab initio 
quantum mechanical approaches. 

This interim report describes our current research effort designed 
to continue the development of a reliable, more accurate, and more efficient 
quantum mechanical method necessary for the calculation of the properties 
of clusters of metal atoms. The method deals with the important simplifica- 
tion in the quantum mechanical calculation of metal atom clusters of the 
removal of the metal atom core electrons from the calculations. The method 
being pursued is based on the Effective Core Potential (ECP) approach to 
simulate the effect of the core electrons on the valence electrons. ^ So 
far in this project we have concentrated on the computer programs that will 
allow us to carry out the quantum mechanical calculations, on possible 
simpl ications in these calculations, on the reformulation of the ECP integral 
computation as efficiency is becoming of the essence for these metal atom 
cluster calculations, and, finally, on the continually expanding need to 
refine the ECP's themselves. In the following section we expand on the 
progress made during these last six months, while the Appendix contains 
a reprint of a publication derived from work on this NASA Grant. 
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SCIENTIFIC PROGRESS 


Evaluation of Matrix Elements of the ECP 


The interface of the program to evaluate the SPDF ECP molecular 
integrals with Komornicki's programs has proven to be straightforward. A 
comparison of the interface codes written independently by both S. Kato and 
ourselves shows essential agreement. The evaluation of the gradients of 
the ECP integrals, in turn, reduces to a matter of calling for special com- 
binations of ECP integrals obtained from the programs thus far interfaced. 
Indeed, S. Kato has recently added this ECP capability to Komornicki's pro- 
grams. One of the advantages of this implementation of the gradient of the 
ECP's over that of the comparable implementation of the SPD ECP's in the 
well-krtown GAUSSIAN 80 package of programs is that the calculation of the 
gradient is not restricted to those cases for which only s and p functions 
appear in the basis set. Plans have been made to use this SPDF ECP en- 
hanced version of Komornicki's programs to investigate the structure of Cs 
atom clusters in collaboration with B. Laskowsky because it has become 
clear that there is a strong common interest in these problems. 

As part of an effort to contribute to the gradient evaluation 

technology, a paper has been submitted and accepted for publication that 

discusses the relationships among derivatives of the molecular integrals 

arising from the transformation properties of the molecular integrals under 

infinitesimal rotations. The paper has been recently published in; J. 

Chem. Phys. 3962 (1981). A preprint of this paper is found in the 

Appendix to this interim report. While the study in this paper finds 

additional simplifying relationships among the derivatives of the integrals, 

analogous to the ones arising from translational invariance as discussed by 
(2) (3) 

Komornicki et ar ' and Dupuis et al ' , the numerical advantages that can 

be derived in practice by use of the additional relationships have not yet 
been explored. 

While the Kahn SPDF ECP integral program was the first program 
with the capability to allow up to f-projection operators in the ECP and 
up to f-type GTO's, it was written with reliability and accuracy in mind 
more than efficiency. A recent independent version of the ECP integral 
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program, written by McMurchie and Davidson, is reported to be slower by a fac- 
tor of 2 . Thus, the Kahn SPDF ECP program remains the most efficient one of 
its kind despite the absence of the deliberate pursuit of efficiency ir its 
design. The major source of inefficiency arises from the neglect of the use 
of the shell structure of the GTO's. This causes the computation of ECP 
integrals sometimes to be of comparable length to that of the two-electron 
integrals (evaluated over shells). This, however, has been of no consequence 
i:i past applications of the ECP method because the basis sets in use with the 
F.CP method were so reduced over standard basis sets that the overall computa- 
tion times had become trivial to begin with. Nonetheless, as the size of the 
molecular problems studied with the ECP method becomes larger, as in the 
study of metal atom clusters here considered, the basis sets again become 
large. The need to improve the efficiency of the SPDF ECP integral program 
for these applications by reorganizing the calculations around the shell ap- 
proach is therefore increasingly imperative. 

The organization of the SPDF ECP integral calculation around the 
shell approach is found not to be straightforward. The general ECP integral 
formulae have been shown by Kahn et al^ ^ to be as follows: 
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The obvious problem is that these formulae do not exhibit the needed dis- 
tinction of terms as being common to the components of a shell from those 
that change from component to component. 

The first step taken to reexpress the integral formulae in a 
more shel 1 -oriented structure was to find an alternative indexing system 
for the basic Gaussian-type orbitals (GTO's). Conventionally, each GTO is 
indexed by the set of "quantum numbers" ( ‘‘'j . ) that determine the 

angular dependence of the GTO as in 

. 

An alternative indexing system is needed that is more explicit about the 
shell characteristics. It is found in the new set of "quantum numbers" 

( "U , X t ji ); these are simply related to the old t IV ) 

"quantum numbers" by ^ 


tv a ifvt, fiVj + (4a) 

Yl«, + (4b) 

The first index, , is but the shell angular-momentum quantum number, 
and it remains the same for all components of a shell. One can range over 
all the "quantum numbers" ( ) of a given shell just by changing 

the % and u. "quantum numbers" as follows 
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The simplicity and order in the ( U , X , ^ ) indexing system over that 
of the ( ^ ^ ) indexing system becomes apparent upon even a casual 

inspection of this f-shell case. The simplicity of the new indexing system 
is demonstrated, for example, in the ease with which the number of components 
of a shell can be calculated, i.e. 
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ihe properties of the GTO's relevant to the ECP integrals must 
now be found in terms of the ( tl- , 1 , ^ ) index system to enable the 
reformulation of the integral formulae. Consider in particular the angular 
dependence factor 
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(8) 

The translation of , 

*c = 

or J from center C onto center A is 
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The expansion of the angular dependence factor ^ from center C onto 
center A induced by this translation can be worked out straightfornvardly 
for any particular case given ( U , 1 , u ). For example, for (2, 2, 0 ) 
one finds 


The general expression, however, is considerably more difficult to find. 

We have found a general expression that shows the expansion terms in the 
new ) "quantum number" system, as follows 

Xciti-Z £ £’ X 01) 
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a -ywMt (Oj 4vf l+tt) , (12a) 

X, = (^,^1 . (12b) 
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In this formula, the coefficient of the expansion is just a product of 
binomial coefficients, 
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The integral formulae may now be reexpressed in terms of the new 
GTO ( H , X . It ) indexing system. Let us first define the basic radial 
integral s 
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The ECP matrix element of the first kind , apart from GTO normalization 
factors, is 
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where tie c'’tffic1ent of the I radial integral is 
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The ECP matrix element of the second kind , apart from GTO normal ization 
factors, is 
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where the coefficient of the J radial integral is 

C**' » E 

r.t;r r 

'{'-M 2 


i »' 

n'+i 


( 21 ») 


(21b) 



9 


I't (t) 
y n-h' 


f *n 


^+1 < ^ i tX 


(21c 


and tne basic 


C. .1 gi «« 


elements are defined as 


(t) 

t:.,.^--z I 6^,. xi«?> . 

i! ft 


I r 6 7 X (*-»’ 
iri; ftV ' 

I LW r l/W) i 1 1, ,X.- i, AA j? .X.,.,., ? Att 

fc'sC i,w' ■!%(“ {>»• -mst SL ^ A 


(22) 


The are the real spherical harmonics. 

The new ECP integral formulae, Equations (17) and (20), replace 
the old formulae given by Equations (1) and (2). The principal result in 
Equations (17) and (20) is that the I and J radial integrals are shown to 
be common to all the components of a shell. Because radial integrals are 
often quite time consuming to calculate, this result alone shows these new 
formulae to yield the desired simplification that is the promise of the 
shell approach. All the terms that depend on the detailed geometry are 
collected in the C coefficients; these coefficients appear to be fully depen- 
dent on the properties of the individual components, although we have not 
yet concluded the analysis. 

One very important simplification has come to light upon exten- 
sively exercising a computer-driven loop through the various angular integrals 
appearing in both the E and E elements. Equations (19) and (22): the majority 

C-. these angular integrals are zero! The remaining task is to identify pre- 
cr:.. all the selection rules for these angular integrals so that only the 
non-vanishing terms in the D and E elements are ultimately evaluated. After 
the analysis is complete, it is planned to implement these new formulae. This 
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will require the complete revision of many of the routines in the current 
SPDF ECP integral program since the structure of the formulae is completely 
different from the previous formulation. The basic interface with Komornicki's 
program, however, is expected to remain essentially unaltered, thereby making 
this upgrade compatible with the work already done. 

Improvements in the ECP's 

A new unified approach to the calculation of the "shape-and 
Hamiltonian-consistent" ECP's has been derived; it contains the "shape- 
consistent" approach of Christiansen et al.' ' as a special case, and it 
deals directly with the essential orbital properties that in turn ensure the 
matching of the valence-valence interactions in the approach of Rappe et al.^^^. 
The foundation of this unified approach lies in the formulation of the atomic 
energy in terms of moment-accumulation functions as presented by Kahn^®^ at 
the recent NRCC-sponsored workshop or Effective Potential Methods. More 
details will be presented in the final report on this project. 
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Relationships among derivatives of the integrals in the calculation 
of the gradient of the electronic energy with respect 
to the nuclear coordinates 

Luis R. Kahn 

Physical Metallurgy Section 
Battel le, Columbus Laboratories 
Columbus, Ohio 43201 
(Received 


ABSTRACT 


The use of the translational invariance property of molecular integrals 
is extended by its combination with the closely related transformation 
properties under infinitesimal rotations. This further reduces the number 
of truly linearly independent derivatives of the molecular integrals that 
require explicit calculation in obtaining the gradient of the electronic 
energy. The detailed linear equations that yield the linearly dependent 
derivatives in terms of the "independent" derivatives are presented. 
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I. INTRODUCTION 

The usefulness of the gradient of the molecular electronic energy 
with respect to the nuclear coordinates ranks close in importance to that 
of the electronic energy itself. The gradient yields an immediate sense 
of the local topography about any point on the potential energy surface 
for a chemical reaction. It enables, moreover, the efficient and systematic 
convergence to the points at which the energy surface is stationary with 
respect to variation of the nuclear coordinates. This is especially 
useful in multidimensional cases in which physical intuition may be of limited 
value in helping to locate these critical points. 

The techniques for the accurate calculation of the gradient of the 
molecular electronic energy have been actively developed and refined in 
recent years. The computational effort required to calculate the 
gradient, notwithstanding these recent important advances, is not a minor 
one. Both Komornicki et al^^^^ and Dupuis et al^^^^, for example, find 
that, on the average, the calculation of the gradient requires about 2-1/2 
times the computational effort of the calculation of the energy itself. 

Apart from simple overlap Integrals, the expression for the gradient 
of the electronic energy differs from that for the energy by the replace- 
ment of the integrals over atomic basis functions by the derivatives 
of these integrals with respect to the nuclear coordinates. Consider, 
for illustrative purposes, the case of a closed-shell system in the single 
determinant approximation. The electronic energy is^^*^^ 


( 1 ) 
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The c. vector contains the expansion coefficients of the i'th molecular 
'ui 

orbital (MO) in the atomic basis functions. The summation in Eq. (1) 
ranges over all occupied MO's. The elements of the matrix are the 
one-electron integrals, and is the matrix of Coulomb and exchange 
two-electron integrals in the well-known form 

(2) 

Sf 

where the elements of the supermatrix B are 

®y,v;»c,X " J [w^Ivk] (3) 

The Tuvlicx] term, for example, is a basic two-electron integral over 
atomic basis functions. Let the a'th component (a = x,y,z) of the 
coordinates of the k‘th atom be denoted as P. . The partial derivative 

K f Ol 

of the electronic energy with respect to P. has been shown by Moccia^^^ 

K ,a 

(8) 

and Fletcher' ^ to be given by 



( 1 )' 




( 4 ) 


where 


and 


(1) 


( 0 ) 


k,a 
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i 




(5) 

( 6 ) 


and 




(7) 


The elements of the S matrix are just the overlap integrals. Thus, 
the calculation of the gradient of the electronic energy is seen from 
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this example, to revolve about the calculation and assembly of the derivatives 
of the one-electron integrals, the derivatives of the two-electron integrals, 
and the derivatives of the overlap integrals, Eqs. (5), (6) and (7) re- 
spectively. 

In this paper we discuss relationships that exist among the various 

derivatives of each integral required in the calculation of the gradient 

of the electronic energy. These relationships derive from the spatial- 

transformation properties of the integrals independently of the existence 

of spatial symmetry in the molecular system as a whole. The specific 

spatial-transformation properties of interest here derive from the molec- 

( 22 ) ( 22 ) 

ular integrals under infinitesimal translations^ ' and rotations,' ' The 
use of the property of translational invariance alone was first proposed 
by Komornicki et al^^^ and later extended by Dupuis et al.^^^^ Whereas the 
translational invariance property has been used to obtain significant simp- 
lifications in the calculation of the gradient of the electronic energy, 
no use has been made of the closely related transformation properties under 
infinitesimal rotations, Dupuis et al have, in passing, commented on the 

special condition of rotational invariance when only s-type basis functions 

fl5l 

appear in the integrals but made no use of it.' ' We show in this work 
the additional relationships among the derivatives of the integrals that 
can be derived from the combined use of the properties of the integrals 
under infinitesimal translational and rotational transformations. We find 
that the six conditions that result from the transformation properties 
taken together reduce the number of linearly independent derivatives from 12 
to 6 in the case of a four-center integral, and from 9 to 3 in the case of a 
three-center integral. In a two-center case, the number of linearly indepen- 
dent derivatives is reduced from 6 to 1. The nu»»^er of derivative integrals 
that require explicit calculation in the evaluation of the gradient of the 
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electronic energy are found to be considerably reduced, in principle, by 
the linear relationships among them. These linear relationships are found 
from the transformation properties of molecular integrals under infinites- 
imal translations and rotations. 


II. RELATIONSHIPS AMONG THE DERIVATIVES OF A MOLECULAR INTEGRAL 
WITH RESP~ECT TO THE NUCLTEAR COORDIN^ 


Let us denote an integral by the symbol I and, although we are concerned 
fit most with a four-center case, let us assume, for the sake of generality, 
that the integral is an N-center case. We indicate this by writing 
^^^1* ^2’***^N^* integral must be translationally invariant. There- 
fore, we may write 




^2 .vn' ' 




^ 2 - 


T P 


-P 


( 8 ) 


where T is an arbitrary three dimensional vector. Let T be an infinitesimal 
translation vector, and expand the right side of Eq. (8) in a Taylor 
series to obtain 


^ 


^ ^Jl» ?2 ’****Sn^ 


E T 

a=x,y ,2 


N 

E ^ 


(9) 


>a 


Equating linearly independent terms results in conditions on the partial 
derivatives of I to all orders. The condition on the first-order partial 
derivatives of I, the same as previously discussed by Komornicki et al 
is found to be 


N 

y - 

k=l ^^k,a 


0 


(10) 


a * x,y,z 



A- 6 


We consider next the additional conditions on the first-order partial 
derivatives of I resulting from the transformation properties of I under 

A 

Infinitesimal rotations. We define the infinitesimal rotation, R, by 


RP 


k,a 


^a.e) 


e 


x,y,2 


( 11 ) 


where e 


a, 6 


e“x,y,z 

are the infinitesimal matrix elements of an antisymmetric 


- e. 


a, 6 6»a 

The integrand of the I integrals depends on the nuclear coordinates 
{P|^} through both operators and basis functions, each centered on one of 
these nuclear coordinates. The nuclear attraction operators are the 
only operators explicitly dependent on the nuclear coordinates in the 
integrals appearing in the calculation of the electronic energy. Since 
these operators are spherically syninetric about their respective nuclei, 
the transformation properties of the I integrals upon an infinitesmal 
rotation are solely determined by the transformation properties of the 
basis functions. 

We take the basis functions to be Cartesian Gaussian basis functions 

\.n .n 
X jf z 

where N{;,r^/iy,n^) is the normalization constant. These functions are 

usefully grouped by shells. The functions in the same shell have the 

property that, in addition to the exponential parameter Ci the sum 

n • n +n +n is the same. The transformation properties of these basis 
X y z 

functions upon an infinitesimal rotation, (11), are 

'o.e * 


( 12 ) 


( 25 ) 


( 13 ) 


X’ y' z 


• • • 


( 14 ) 
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where 


6G 
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m t 
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A G 

x,y 

(15a) 
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A G 
x ,2 

(15b) 

6G 

2 ,y 
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\.2®nj^.ny-l,n^+l 

(15c) 
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2n„ + 



^a,6 

•"»( 

2n - l) 

a.B * x,y,z 

(15d) 


Equations (14) and (15) show that, upon an infinitesimal rotation, a basis 
function transforms into other functions of the same shell. The integrand 
of I consists, apart from the usual operators, of a product of these basi:< 
functions. The transformation properties of the integral I under an in- 
finitesimal rotation are, consequently, found to be 


*zz 


o 4^ •Pm)"’’*** 

a»6 a,$ 'ol '^2 ^<N 


(16): 


where 



(0 

dl 


(17) 


The summation in Eq. (17) ranges over all basis functions appearing in 

the integrand of I. The temi is defined as an I integral in which 

the I'th basis function, , has been replaced by a function fiG 

"x'"y"2 (.) 

defined by Eqs. (15a) - (15c). Consequently, each Integral 6V ' is just 

01 f P 

a simple linear combination of two I-type of integrals, as follows 


*t(t) . A j(l) . A 

y.x y,x^nj^*»^l,n^-l ,n^ x,y*nj^-l,ny+l,n^ 

xl(t),« .A !(*•) 

2.x ^'z,x 


(18a) 


(18b) 
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®^2»y i.y ^ytZ^n^.n^-l .n^+l 

/.\ fp) 

The notation I'y' , analogously to its use in 6 l denotes 

n jj»n ^,n ^ a,e 

an I integral in which the I'th basis function, , has been 

%»"y*"z 

replaced by another basis function, G., . Note that the re- 

n ^.n ^,n ^ 

placement functions belong to the same shell as the original function. 
The effect of an infinitesimal rotation can also be expressed as 


(18c) 


•IRPr “Sn> • '<?i * ^2 * t?2 Sri * *i> 

Expanding the right side of Eq. (20) in a Taylor series, and using the 
antisyrmetry property of the ; matrix, one obtains 


R]^2 


♦IE 


e 


a, 6 



k,6 


- P 


k • Cl 


k»a 



( 20 ) 


Wc equate each linearly independent set of terms 1 n this result with 
those obtained in Eq. (16). This yields the conditions on the first- 
order partial derivatives of I that derive from the transformation 
properties under infinitesimal rotations. One obtains 



> 4 

k,o JP, 


6l 


k,e 


0,6 


0,6 ■ ytx; 2 ,x; 2 ,y 


(21) 


It follows from Eqs. (18a) - (18c) that the right side of Eq. ( 21 ) 

Is zero in the special case when all the basis functions appearing 

In I are s-type functions, n -n -n ■ 0. Equation (21) reduces in 

X y 2 

(15 ) 

that special case to a condition of rotational invariarice^ . 
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Taken together, Eqs. (10) and (21) constitute the six conditions relating, 
in general, the first-order partial derivatives of an arbitrary molecular 
integral I to each other. There are many possible ways to combine Eqs. 
(10) and (21). One particular combination to be used later Is 



.B 




■ 6l 


0,6 


0,8 “ y.x;z,x;z,y (22) 

The gradient of I has 3N components. One can always obtain 
(3N-6) of these by independent explicit calculation and then, it would 
seem at first, that one could obtain the remaining six components by 
straightforwardly solving the system of linear equations. A closer 
inspection reveals, however, that the rank of the coefficient matrix of 
the reduced system of linear equations depends on the choice of the (3N-6) 
"independent" components. We therefore discuss three imoortant applications 
of Eqs. (10) and (21) separately. 


A. Gradients of Two-Center Integrals 

The gradient of the two-center Integral has six components. 

Since there are also 6 conditions, Eqs. (10) and (21), relating the six 
components to each other, one would at first conclude that these equations 
alone fully define all the components of the gradient. A reduction of 
the coefficient matrix of the sysiem of linear equations to canonical form, 
however, shows that the rank of this matrix^^^^ Is actually only 5. Thus, 
one of the six components of the gradient may be considered truly inde- 
pendent. We choose this one to be »I/aP, for example, since the choice 

I fX 



Is inwaterial in this case. Assu»n1ng that ^ the remaining 

I f ^ t X 

5 components of the gradient are obtained as 







il 


♦ fil 


1.x 



o ■ y,2 




a ■ x.y.z 


B, Gradients of Three-Center Integrals 

The gradient of the three-center Integral has nine components. 

The rank of the coefficient matrix of the linear equations relating 
these nine components can be shown to be 6. We exclude for the purpose 
of this discussion the special case of col linear geometries. In this 
special case the rank of the coefficient matrix is 5. Moreover, 
because of the special geometry, the number of implicit independent 
variables is reduced as well. Excluding the collinear case, three of 
the nine components of the gradient may hence be taken as "independent" 
components. Not every choice of these three components, however, yields 
a system of linear equations of rank 6. Examples of choices of the 
three components that lead to a rank lower than 6 are 



and 



k- 1.2,3 (25) 


a • x,y,z (26) 


A choice of three Independent components that preserves the rank of 6 
should include at least two different centers, and, if three different 
centers are present, at least two different types of cartesian 
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components. Of the many possible choices that satisfy these require- 
ments, we shall use the following 



(27) 


The equations that follow are therefore only representative of 
the linear equations determining a set of "dependent” gradient 
components in terms of a counterpart choice of "independent" 
components. 

It is convenient at this point to define a notation that will 
allow us to express in compact form the remaining equations. First 
we define the elements of a 3 x 3 matrix A, „ as 


■ ^N,z' 




-(^.z - ^N.z’ 


- ’’n.x) 


^’’z.y ■ ^N,y*\ 

(^.x - ’’n.x) ] ( 28 ) 


where N is equal to 3 in the three-center irtegral case. Second, 
we define the components of a three dimensional vector ^ as 


( 


3l 3l 3l \ 

^ 


(29) 


Note, for example, that Expression (27) gives the components of a 
9. , vector, and that 9. , + 9« , = 9, , + g„ «. Lastly, we define 

'v.l.Z 'V/l.Z -x/Z,! '\xl,l 'b2,2 '' 

the components of another three dimensional vector 6 as 





( 30 ) 
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The elements of this vector have been in turn defined in Eqs. (17) and 
(18). 


The ''independent'' components of the gradient have been chosen 
as the three elements of the g, , vector, Expression (27). The re- 

<\,i 

maining six components are contained in the vectors ^2 -\ 93 3* 

These vectors are obtained in turn as the solutions to three simultaneous 
linear equations. Assuming P9 , ^ Po these equations, in matrix 
notation, ere 

'^ 2,1 ? 2.1 ‘ I 

?3.3 ' ‘ ?1,2 " §2,1 

where the vector on the right side of Eq. (32), f, is defined as 

tl,2 §1,2 i 


f = 


(32) 

(33) 

(34) 


C. Gradients of Four-Center Integrals 

The gradient of the four-center integral has 12 components. 

The rank of the coefficient matrix of the linear equations relating 
these 12 components can again be shown to be 6. Therefore, six out 
of the 12 components of the gradient may be taken as the "independent" 
components. A convenient (but not unique) choice of "independent" 
components that preserves the rank of 6 in the reduced linear equations 
are the elements of the two vectors -j and ^ y Assuming P3 ^ ^ ^4 z 
and excluding again the special case of col linear geometries, the 
equations giving the remaining components of the gradient are 

^3.2 ?3,2 ' I 


(35) 
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?4.4 ” ‘ ?1.1 " ?2.3 ’ ?3,2 

The vector ^ on the right side of Eq. (35) is defined as 

I ' • ^ 2 , 3 ? 2,3 ‘ ^ 1 . 1?!. 1 * ^ 


(36) 


(37) 


The various ^ matrices in this case contain N set equal to 4. 

In the general case, the “independent" components may be taken as 
the vectors ^_2 f^_<| and k * 1 , N-3}. The remaining com- 
ponents of the gradient are contained in the vectors l ^_2 

and ^ 1 ^. Assuming z ^ z* components are obtained from 

the equations 

^N-1 ,N-2 ?N-1 ,N-2 ' ^ (38) 


N-3 

^N,N " ■ Sk,k ■ ?N-2,N-1 ' ?N-l,N-2 
where the vector f on the right side of Eq. (38) is 


(39) 


N-3 

J ^ ■ ^-2,N-1 ^N-2,N-1 ■ Z ^k.k ?k,k ”*■ i 

K=1 


(40) 


III. CONCLUSION 

The calculation of the gradient of the electronic energy requires 
the calculation of the derivatives with respect to the nuclear coordinates 
of the integrals over atomic basis functions. The number of these de- 
rivative integrals to be calculated is very large, particularly because 
of the two-electron type of integrals, of which there are more than 
M^/8 in the case of M basis functions. We have shown that the condition 
of translational invariance of the integrals can be combined with the 
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conditions expressing the transfon.iation properties of integrals under 
infinitesimal rotations to obtain equations that yield the “dependent" 
derivative integrals in terms of the "independent" ones. These equations 
are shown to be at most a set of three simultaneous linear equations. 

The existence of these interrelationships, in principle, reduce con- 
siderably the number of derivative integrals that require explicit cal- 
culation. The number of linearly independent derivatives is reduced from 
12 to 6 for a four-center integral, from 9 to 3 for a three-center integral, 
and from 6 to 1 for a two-center integral independently of the existence 
of spatial symmetry in the molecular system as a whole. 
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